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LONG-TERM  GOAL 

Our  long-term  objective  is  to  understand  how  deep  convection,  induced  by  strong  buoyancy 
forcing  at  the  ocean  surface,  influences  the  ocean  circulation  through  convective  plumes  and 
geostrophic  eddies. 

SCIENTIFIC  OBJECTIVES 

Issues  specific  to  oceanic  deep  convection  are  the  relatively  strong  role  of  rotation,  such  that 
convective  Rossby  numbers  appropriate  to  the  vertical  motions  may  fall  below  unity;  the  highly 
intermittent  nature  of  the  forcing,  with  the  strong  impulses  coming  during  the  passage  of 
atmospheric  weather  systems;  the  possible  highly  localized  occurrence  of  deep  convection;  and 
the  subsequent  spatial  redistribution  and  mixing  of  convected  water  by  geostrophic  mesoscale 
dynamics,  which  arise  in  part  from  instabihty  of  the  localized  convection  regions.  Our  particular 
research  projects  are  computational  investigations  intended  to  (1)  quantify  the  extent  to  which 
penetrative  mixing  occurs  at  the  base  of  a  rotation  convective  layer,  sharpening  the  pycnocline 
and  entraining  denser  fluid  from  below;  (2)  determine  how  space/time  intermittence  in  either  the 
surface  buoyancy  flux  or  the  pre-existing  ocean  stratification  and  circulation  relate  to  the  degree 
of  localization  of  deep  convection;  and  (3)  explore  the  circumstances  under  which  localized 
converted  regions  are  persistent  after  the  buoyancy  forcing  abates  or  undergo  erosion  through 
subsequent  horizontal  mixing  and  restratification  of  the  gyre  interior. 

APPROACH 

High-resolution  numerical  simulations  using  the  Boussinesq  model  developed  in  cooperation  with 
our  colleagues  Dr.  Keith  Julien  and  Dr.  Joseph  Weme  of  NCAR  are  carried  out  on  grids  of  256  x 
256  X  129  over  domains  of  2km  x  2km  x  1km  to  investigate  plume-scale  dynamics.  Conditionally 
sampled  composite  techniques  are  used  to  isolate  the  plume  structures.  The  Boussinesq  model  is 
also  used  at  coarser  resolution  over  domains  of  50  km  x  50  km  x  2  km  to  examine  interaction 
between  convection  and  simple  examples  of  geostrophic  circulation. 

For  more  complex  inhomogeneities  in  the  initial  conditions,  larger  horizontal  domains  and  longer 
evolution  times  are  required.  Limited  computational  resolution  then  prevents  resolution  down  to 
the  plume  scale,  and  a  parameterization  of  the  convective  mixing  is  applied  in  these  large 
domains  (up  to  200  km  x  200  km  x  4  km).  The  use  of  an  implicit  gravity-wave  Boussinesq  model 
developed  by  Alistair  Adcroft  at  MIT  enables  the  use  of  longer  timesteps,  while  representing  the 
geostrophic  adjustment  process.  We  are  studying  the  interaction  between  the  vertical 
homogenization  forced  by  the  surface  buoyance  loss  and  the  preexisting  eddy  field,  including  the 
development  of  baroclinic  instabihty,  the  horizontal  homogenization  of  the  mixed  fluid  by  the 
eddy  dynamics,  and  the  possibihty  of  persistence  of  isolated  cores  of  mixed  fluid,  and  comparing 
these  phenomena  with  observations  from  the  field  program. 


1 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number 

1.  REPORT  DATE 

30  SEP  1997  ‘REPORT  TYPE 

3.  DATES  COVERED 

00-00-1997  to  00-00-1997 

4.  TITLE  AND  SUBTITLE 

Deep  Convection  in  the  Ocean 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  California,  Los  Angeles, Institute  of  Geophysics  and 

Planetary  Physics, Los  Angeles, CA, 90095 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a  REPORT  b  ABSTRACT  c  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Office  of  Naval  Research  Space  and  Remote  Sensing 

1997  Annual  Report 

WORK  COMPLETED 

The  following  works  have  been  completed: 

1.  We  have  completed  a  series  of  numerical  simulations  of  convection  forced  by  homogeneous 
cooling  into  a  volume  of  fluid  containing  a  single  geostrophic,  cold-core  eddy,  providing 
extensive  evidence  of  localization  of  convection  by  anomahes  in  the  preexisting  ocean 
stratification.  We  have  examined  the  dependence  of  both  the  degree  of  locahzation  of  the 
convention  and  the  subsequent  barochnic  instabihty  on  the  strengths  of  the  initial  eddy 
circulation  and  stratification  anomaly  and  on  the  size  of  the  eddy  relative  to  the  deformation 
radius.  We  have  performed  a  detailed  analysis  of  the  results,  examining  in  particular  the 
relative  importance  of  buoyancy  fluxes  due  to  a  secondary  circulation  compared  to  those  by 
small  scale  plumes  and  eddies.  These  results  are  described  in  the  article  “Locahzation  of  deep 
convection  by  a  mesoscale  eddy”  (Legg,  McWilhams,  and  Gao)  in  press  in  the  Journal  of 
Physical  Oceanography; 

2.  We  have  completed  an  analysis  of  the  plume  structures  of  turbulent  convection.  Plume 
structures  are  extracted  from  the  turbulent  fields  associated  with  convection,  by  conditionaUy 
samphng  for  maxima  in  the  downward  vertical  velocity,  and  calculating  a  composite  plume. 
Using  the  composite  to  identify  typical  plumes  we  have  obtained  a  census  of  the  plume 
population  and  examined  statistics  including  the  amphtude,  vertical  tilt,  horizontal  radius, 
axisymmetry,  and  vertical  length.  The  composite  plume  is  compared  with  simple  plume 
models  to  identify  the  predominant  physical  balances  including,  in  particular,  the  relative 
importance  of  entrainment  and  detrainment  of  plume  fluid  and  forces  acting  to  enlarge  and 
accelerate  plumes  crossing  the  convecting  layer.  The  results  of  this  plume  analysis  are 
described  in  the  article  “Plumes  in  rotating  convection:  Part  I  Ensemble  statistics  and 
balances,”  Juhen,  Legg,  McWilhams,  and  Weme,  submitted  to  the  Journal  of  Fluid 
Mechanics; 

3.  We  have  developed  several  versions  of  a  lagrangian  particle  code,  to  be  used  in  investigating 
the  response  of  massless  particles  to  convective  flow  fields.  One  version,  in  which  particles 
are  constrained  to  particular  vertical  levels,  similar  to  isobaric  floats,  has  been  used  to 
investigate  the  response  of  floats  to  convection  in  the  presence  of  mesoscale  eddies.  Another 
version,  in  which  particles  are  advected  by  the  fuh  3-D  flow  field  has  been  used  to  investigate 
the  evolution  of  fluid  parcels  in  individual  convective  plumes,  to  complement  the  composite 
analysis  described  above.  These  latter  results  are  described  in  the  article  “Plumes  in  rotating 
convection:  Part  II  individual  plume  hfecycles”  Juhen,  Legg,  McWilhams,  and  Weme, 
currently  in  progress; 

4.  We  have  cohaborated  with  Ken  Fischer  et  al.  at  the  Environmental  Research  Institute  of 
Michigan  by  sharing  surface  fields  from  our  locahzed  convection  simulations,  which  they 
have  used  to  simulate  Synthetic  Aperture  Radar  (SAR)  images  of  the  convecting  fields. 
These  simulations  investigate  the  feasibility  of  using  SAR  to  observe  the  structures  of  open- 
ocean  deep  convection.  An  article  describing  the  results,  and  the  feasibihty  assessment, 
“Modeled  radar  surface  signature  of  deep  ocean  convection,”  K.W.  Fischer,  S.  Legg,  W.H. 
Munk,  R.M.  Shuchman,  R.W.  Garwood,  and  J.P.  Palshook,  has  been  submitted  to  IEEE 
Transactions  on  Geoscience  and  Remote  Sensing; 
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5.  We  have  carried  out  a  series  of  calculations  of  convection  in  the  presence  of  mesoscale 
eddies,  where  the  eddies  are  associated  with  compensating  anomalies  in  Temperature  and 
Salinity,  in  addition  to  the  density  anomaly  previously  studied.  As  the  flow  evolves  in 
response  to  surface  buoyancy  loss,  considerable  structure  on  a  variety  of  scales  is  generated 
in  the  temperature  and  sahnity  fields;  and 

6.  Using  the  imphcit  gravity-wave  model,  we  are  investigating  a  variety  of  convection  scenarios 
in  which  the  initial  conditions  contain  a  field  of  mesoscale  eddies,  with  resultant  interactions 
between  the  eddies  determining  the  eventual  degree  of  homogenization  of  the  flow. 

RESULTS 

The  results  are  hsted  as  the  following: 

1.  Localization  of  convection: 

Our  simulations  of  convection  in  the  presence  of  geostropic  eddies  indicate  that  such 
stratification  anomalies  may  provide  a  mechanism  for  the  long-hypothesized  preconditioning 
of  the  ocean  for  locally  deeper  convection.  Such  anomalies  also  provide  a  mechanism  for 
restratification  of  the  convectively  mixed  region,  as  vertical  convection  is  replaced  first  by 
slant-wise  convection,  and  then  barochnic  instability  with  an  associated  secondary 
circulation.  Instability  eddies  carry  the  dense  fluid  outwards,  causing  the  region  of  deepest 
vertical  buoyancy  flux  to  migrate  outwards.  The  horizontally  averaged  density  field  which 
results  has  hghter  surface  layers,  and  denser  layers  than  in  that  generated  by  horizontally 
homogeneous  convection; 

2.  Plume  structures  at  low  Rossby  number: 

The  plume  structure  identification  algorithm  allows  us  to  examine  the  hfe  cycle  of  a  plume, 
from  emergence  from  the  boundary  layer,  through  its  migration  across  the  convective  layer. 
We  compared  the  properties  of  the  evolution  of  the  typical  plume  vertical  momentum, 
volume,  and  heat  content,  with  an  idealized  entraining  plume  model  such  as  those 
summarized  in  Turner,  1986.  Results  indicate  that  initially  entrainment  leads  to  expansion  of 
the  plume,  but  this  entrainment  quickly  ceases,  and  is  replaced  by  strong  mixing  between  the 
plume  and  environment  which  modifies  plume  properties  with  no  change  in  volume  (see 
Figure  1).  Examination  of  individual  plume  hfecycles  and  the  fate  of  particles  contained 
within  such  plumes  indicates  that  mixing  events  are  associated  with  plume -plume 
interactions,  such  as  mergers,  filamentation  and  shearing,  caused  by  the  cyclonic  vortices 
associated  with  the  plumes. 

Comparison  of  terms  in  the  vertical  momentum  equation  indicate  that  buoyancy  force  acting 
on  the  plume  is  largely  balanced  by  a  combination  of  pressure  drag  and  momentum  exchange 
with  the  environment.  The  resultant  acceleration  is  therefore  significantly  reduced.  Our 
solutions  of  strongly  rotating  penetrative  convection  (Julien,  et  al.,  1996)  have  previously 
indicated  significantly  less  penetrative  convection  as  the  Rossby  number  is  decreased:  both 
increased  detrainment  of  plume  fluid,  and  suppressed  plume  acceleration  are  possible 
candidates  for  this  decreased  penetration; 

3 .  Simulated  Synthetic  Aperture  Radar  ( SAR )  fields: 

The  simulated  SAR  images  reveal  clearly  the  strong  shear  and  strain  in  the  surface  velocity 
fields  of  localized  convection,  associated  with  the  generation  of  cyclonic  vorticity  in  the 
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regions  converging  into  plume  and  barocUnic  eddy  downwelling  regions.  Eddies,  fronts  and 
plumes  can  all  be  identified  and  differentiated  in  these  regions; 

4.  Behavior  of  floats  in  conv  acting  flows: 

Our  investigations  of  isobaric  floats  in  the  presence  of  convection  locahzed  by  mesoscale 
eddies  reveal  a  tendency  for  floats  to  congregate  in  the  convergent  regions  associated  with 
downwelhng  fronts.  These  fronts  are  associated  with  cooler  than  average  fluid  -  hence  the 
ensemble  mean  temperature  measured  by  the  floats  diverges  from  the  Eulerian  mean,  an 
important  consideration  when  interpreting  observations  made  by  such  isobaric  floats.  In  flows 
with  considerable  eddy  structure,  the  dispersion  of  floats  can  also  reveal  the  process  of  lateral 
exchange  and  restratification;  and 

5.  Temperature  and  salinity  fine-structure: 

The  addition  of  compensating  Temperature  and  Salinity  signatures  to  the  eddy  field  provides 
a  means  of  distinguishing  water-masses  originating  from  the  surface,  the  eddy  interior  and  the 
ambient  flow.  Within  a  region  of  approximately  uniform  density,  fluid  having  its  origins  in 
these  different  areas  may  be  found,  leading  to  significant  small-scale  structure  in  the 
temperature  and  sahnity  fields,  as  a  result  of  the  incomplete  nature  of  mixing  (see  Eigure  2). 
Eor  example,  within  an  eddy  generated  through  instability  of  a  weakly  stratified 
preconditioning  eddy,  fluid  of  a  particular  density  consists  of  a  combination  of  water  from  the 
preconditioned  interior  uninfluenced  by  the  surface  forcing,  water  originally  from  the  surface 
made  denser  through  the  surface  cooling  and  advected  downward  along  a  narrow  front,  and 
water  from  the  deeper  layers  of  the  exterior  fluid,  upweUed  in  broader  regions.  Geostrophic 
eddy  activity  stirs  these  different  water  masses  together,  generating  interleaving  filaments  of 
each  water  type.  The  horizontal  scales  of  variabihty  are  being  compared  with  those  noted  in 
the  observations  (Pickart,  private  communication),  and  timescales  for  persistence  of  these 
structures  estimated. 

IMPACT/APPLICATION 

Our  research  on  convective  dynamics  should  improve  understanding  of  both  the  observed 
features  of  convection,  and  the  net  effects  of  these  small-scale  features  on  water  mass 
transformation.  This  understanding  provides  a  firmer  basis  for  predictions  of  sub-polar 
circulations  and  tracer  distributions  and  the  global  thermohaline  circulation. 

TRANSITIONS 

In  September  1997,  Eegg  moved  to  Woods  Hole  Oceanographic  Institution,  to  take  up  an 
assistant  scientist  appointment.  Our  PYI998  continuation  of  this  grant  will  therefore  be  divided 
between  the  two  institutions  (UCEA  and  WHOI)  with  Eegg  taking  responsibihty  for  the  WHOI 
component.  Eegg  submitted  a  separate  proposal  for  this  component  in  August. 

RELATED  PROJECTS 

The  following  is  a  list  of  related  projects: 

1 .  ONR  Deep  Convection  ART. 
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Oiu'  simulations  of  convection  localized  by  one  or  more  mesoscale  eddies  provide  scenarios 
for  cou^arisou  with  observations,  both  of  the  small  scale  velocity  and  temperatiue  fields,  and 
remote  sensing  signatiues  of  the  chimney  scale  circulation  and  barochnic  instabrhty  eddies. 
These  simirlations  encoitrage  the  evahration  of  the  pre-convection  eddy  field  for  comparison 
with  the  later  localized  convection  regions.  Oin  plirme  stnrctiue  analysis  provides  detailed 
statistics  on  plirme  property  distributions,  for  comparison  with  data  fiom  observations; 

2 .  ONR  Morim  Boundary  Layer  ARl: 

Oiu'  stirdies  of  deep  convection  driven  by  siuface  cooling  overlap  with  stirdies  of  the  oceanic 
and  atmospheric  planetary  boimdary  layers.  We  are  engaged  in  a  long-standing  partnership 
with  Chin-Hoh  Moeng  and  Peter  Sulhvan  of  NCAR,  making  Large  Eddy  Simulations  of 
boimdary-layer  tmbuleuce  and  comparisons  with  the  field  observations  made  in  this  program; 
and 

3.  NSF  High  Performance  Computing  and  Communications  (HPCC)\ 

Oiu'  research  into  the  small-scale  featiues  of  tiubuleut  penetrative  convection  continues  to 
provide  interaction  with  oiu  colleagues  Keith  Juhen  and  Joe  Werne,  as  part  of  the  University 
of  Colorado/NCAR/UCLA  NSF  HPCC  Grand  Challenge  project. 
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Figiue  1.  (a)  The  entrainment  coefficient  a  = —  where  Veis  the  entrainment  velocity,  and  W  is 

W 

Ux 

the  vertical  velocity  of  a  plume,  and  (b)  the  mixing  coefficient  ttx  =  —  where  U*  is  the  mixing 

W 

velocity  and  W  is  the  vertcal  velocity  of  a  plirme. 


Entraimueut  leads  to  an  expansion  of  a  plimie,  while  mixing  modifies  plirme  density 
with  no  volume  change.  Results  are  shown  for  simulations  of  convective  plirme 
ensembles  at  several  different  values  of  convective  Rossby  number.  Note,  the  reduction 
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in  enteraiument  coefficient  for  low  Rossby  number.  Mixing,  acting  so  as  to  cool  hot 
plumes,  remains  significant. 
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Figiue  2.  Root  mean  square  values  of  (a)  aT  +  ps  and  (b)  pS  -  aT  =  (p  -  po)  /  po. 


T  is  teiiqjeratiue,  S  is  salinity,  p  is  density,  po  is  a  reference  density,  a  is  the  theimal 
expansion  coefficient  and  p  =  •  p  /  •  S.  The  ims  is  defined  as  Xm^  =  —  A”)  where 

A'  is  the  asimiithal  average  of  X.  Fields  are  initially  asimuthally  symmetric,  and  are 
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shown  as  a  function  of  radial  and  vertical  distance.  A  dense  asimuthaUy  symmetric  eddy 
associated  with  temperature  and  salinity  anomalies  is  forced  by  a  uniform  surface  heat 
loss.  Slant-wise  interleaving  and  eddying  associated  with  barochnic  instability  lead  to 
strong  compensating  variability  in  T  and  S.  Hence  variability  in  aT  -i-  (3S  is  significantly 
greater  than  variabihty  in  (p  -  po)  /  po.  Parameters  are  shown  non-dimensionalized  by 
the  surface  buoyancy  loss. 
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